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Of 



We present an analysis of selection biases in the M bh — a relation using Monte- 
Carlo simulations including the sphere of influence resolution selection bias and 
a selection bias in the velocity dispersion distribution. We find that the sphere 
of influence selection bias has a significant effect on the measured slope of the 
M bh - a relation, modeled as beta intrinsic = -4.69 + 2.22(3 measured , where the 
measured slope is shallower than the model slope in the parameter range of /3 > 4, 
with larger corrections for steeper model slopes. Therefore, when the sphere 
of influence is used as a criterion to exclude unreliable measurements, it also 
introduces a selection bias that needs to be modeled to restore the intrinsic slope 
of the relation. We find that the selection effect due to the velocity dispersion 
distribution of the sample, which might not follow the overall distribution of 
the population, is not important for slopes of /3 ~ 4-6 of a logarithmically linear 
Mbh — c relation, which could impact some studies that measure low (e.g., (3 < 4) 
slopes. Combining the selection biases in velocity dispersions and the sphere of 
influence cut, we find the uncertainty of the slope is larger than the value without 
modeling these effects, and estimate an intrinsic slope of (3 = 5.28~^q^. 



Subject headings: black hole physics — galaxies: general — galaxies: fundamen- 
tal parameters — methods: statistical 



Introduction 



The relation between super-massive black hole (SMBH) masses, Mbh, and their host 
galaxy velocity dispersion, a, is one of the central themes in galaxy and AGN studies. The 
number of reliable SMBH estimates has more than quadrupled since the early stages of 



investigation (e.g.. IFerrarese fc 
containing ~ 64 galaxies (e.g., 



Merritt 



2000 



Gebhardt et al. 



Graham et al. 



2011 



2000) , with th e latest samples 



McConnell et al.ll201ll ). Even with an 
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i ncrease in sample s ize, the relation remains tight, e.g., with a scatter of only 0.43 dex 
( 1 Graham et al.ll201ll ). tantalizingly suggestive of a strong link between SMBH and galaxy 
formation, and indicative of AGN feedback processes. Many theoretical predictions and 
models are based on this relation, describing t he interaction bet ween SMBH and the bulge of 
the g alaxy through energy-dr iven winds (e.g.. ISilk &: Reeslll998l ) or momentum-driven winds 



Fabianlll999t lKingll2003l ). The M^ — a relation has been used to supp ort many theories 



,e.g. 

conn e cting the central black hole to other host properties of the galaxy fe.g-.lFerrarese fc Ford 



20051: 



Adams et al 



Okamoto et al 



suggested (e.g 



Graham et al 



200 A. where connections as far as with dark matter halos have also been 
Booth fc Schavd 120101) . There are also studies (e.g 



2003; Merloni &: Heinz 2008; Bower et al. 2006; Pi Matteo et al. 



2008 



2011 



Guftekin et al. 2009 



Beifiori et al.ll2012f ) that examine the possibility of the Mbh — a relation 



varyi ng with the morphology of the galaxy and evolution with redshift (e.g. 
2006h . 



Peng et al. 



The slope and normalization of the Mbh — a relation are important parameters for 
constraining the nature of the feedback mechan ism from the blac k hole. Various models 



and simulations generally predict M bh oc <r 5 (e.g., [Silk &; Reeslll998l ) if the feedback is in the 
mechanical form and M m, oc a 4 if the feedback is dominated by m omentum exchanges (e.g., 



Fabianlll999l : lKindl2003l : iGranato et al]l2004 Murray et al.l 120051 ). The normalization of the 
relation constrains the feedback efficiency, where rj ~ 0.05 is generally required for energy 
feedback models and r] ~ 1 for momentum feedback models. Thus, momentum feedback 
models require that the majority of the black hole growth occurs in the obscured phase, 
which can be tested by future hard X-ray surveys. Energy feedback models suffer from 
cooling problems where the feedback energy is lost by Compton scattering with photon s 
from AGN or other radiative cooling mechanisms (e.g., iKingi 120031 ; ISilk fc Nusserl 120101 ). 
Even the 5% feedback efficiency for energy feedback models is difficult to achieve based on 
current observations of AGN winds. Observationally, there are a number of attempts to 
constrain the Mbh ~ c relation (Table 1), where the slope is measured between 4-5 but with 
large enough errors and variety in measurements that it is unclear which feedback model 
drives the relation. The average slope of major studies is (3 = 4.7, which is inclined towards 
an energy-driven feedback scheme, but the samples used to measure the Mbh — & relation 
are still under close scrutiny. 

Considering the importance of the Mbh ~ c relation, it is essential to thoroughly in- 
vestigate the relation. Authors have previously exa mined the effect of the uncertainties in 
black hole and velocity dispersion measurements (e.g.. lNovak et al.ll2006l ; iMerritt fc Ferrarese 
20011). as well as explored selection biases in choosing samples of quasars versus nearby AGN 
(jSchulze fc Wisotzkil 120111 ). The treatment of errors arguab ly can change the measured slope 
of the Mbh — o relation, and as IMerritt fc Ferraresd (120011 ) discuss, the limited sample size 
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can also impact the relation. Recently, questions were raised whether the fitting results of 
the Mbh — a relation are affected by selection effects and even whether the relation itself could 
be an artifact of selection effects. Studies have tended to use resolution of the sphere of in- 
fluence, Tj, around the SMBH as crit eria for inclusion of a galaxy in t heir sa mple as a proxy 



for rel iable mass measurements (e.g.. iFerrarese fc Fordll2005l ; lHull2008l ). while iGultekin et al. 



(120091 ) argue that no such criteri a is needed. Restricting samples to galaxies where the 
sphere of influence, r{= GM bh /a 2 (|Peebleslll972} ). is resolved provides a lower bound to what 
is observable in the Mbh — a plane. Therefore, directly fitting thes e filtered samples ca n 



result in parameters different from those of an intrinsic relation (e.g., IGultekin et al.ll2009l ). 
Furthermore, a relation measured from a sample using t his criteria might even arise sim- 
ply from setting the lower limit for Mbh — a data pairs. iBatcheldorl (120 lOh addresses this 
issue by applying both the rj resolution cut and the observed upper limit plus scatter in the 
Mbh — cr plane and can reproduce the M bh — a relation in simulations. In particular, the 
author uses observed a and simulated Mbh, evenly distributed logarithmically over 10 1 M Q 
to 10 10 M Q , and makes an cut based on an Hubble Space Telescope (HST) resolution of 
O'.'l and uses the upper scatter of the Mbh — cr relation itself to create an upper bound. This 
simple approach provides a physically different interpretation of the observed Mbh — cr rela- 
tion: instead of a correlation, there is an envelope that sets an upp er bound fo r Mbh given 
cr, which may represent a more realistic model for black hole growth (lKing!l2010l ). The lower 
bound therefore only arises from selection effects in this model. However, this model predicts 
a distribution of detection rates o f black hole masses that are arguably not supported bya 
observations (IGultekin et al.ll201ll ). 



Aside from the sphere of influence selection bias, there can be another selection effect 
in the a distribution that is also worth exploring. Ideally, the relation between Mbh an d cr is 
best constrained if the individual parameters' measured distributions in a small sample follow 
the general distributions from a large sample. Although it is impossible to compare to the 
true M b h distribution, we find differences between the observed a distribution in the M h h — cr 
sample and the distribution from a large sample of galaxies. As shown in our simulation, 
this difference can also change the best-fit values for model parameters, especially when the 
underlying relation is non-linear or skewed. In this paper, we model both the effects of 
the sphere of influence resolution cr iteria and selectio n effects caused by a distributions. In 
addition, continuing the argument of lBatcheldorl (120101 ). we extend the question of whether an 
upper bound can be naturally produced for the observed Mbh ~ cr distribution by simulating 
two more general mass distributions of SMBH masses independently of the cr of galaxy bulges. 
Finally, most previous studies use various direct linear regression techniques to measure the 
parameters of the Mbh — cr relation, which make it difficult to model various selection effects, 



and can certainly impact the final measured slope (IMerritt fc Ferraresdl200ll ). In this paper, 
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we use a Bayesian Monte Carlo approach to robustly measure the parameters of the Mbh — g 
relation including selection effects. 

In Section 2, the samples used are described. Section 3 discusses methods and results, 
followed by conclusions in Section 4. Throughout this paper, unless otherwise noted, we 
assume a cosmology of H =70 km s -1 Mpc -1 , Qm=0.27, and Q\=0.73. 



2. The Samples 



2.1. Observed M bh — a Relation Samples 

Although there are several different compilations of measurements for the Mbh — ° 



relation (e.g., iFerrarese Sz Fordl l2005t iGiiltekin et al.l 120091 ; iGraham et al.l 1201 ll ). most of 



them share and/or repeat data on the same galaxies. Table [T] lists the major samples from 
the literature, with values for the parameters in the Mbh — cr relation when written in its 
standard form, \og(Mbh/ M Q ) = a + /31og(<r/200km s _1 ). The variety o f samples means that 



selecting one with which to work is not a simple process. We choose IGraham et al.l (12011 



hereafter GR11) as the main samp le, with 64 galaxies. The sample contains a ll of the galaxies 
from both IFerrarese fc Fordl (120051 . hereafter FF05) and lGiiltekin et al.l (120091 . hereafter G09), 
with the most up-to-date values at the time. GR11 excludes four galaxies present in other 
samples: IC1459, for which gas and stellar dynamical models differ; NGC2748, for which 
dust is an issue; NGC4594, for which the re is no 3-integral m odel; and NGC7457, where the 
AGN/nuclear core distinction is blurred (IGraham et al.ll201ll ). GR11 also includes a further 
16 galaxies not previously included in the Mbh — & literature. 

To keep the sample volume limited to within 100 Mpc, we eliminate four galaxies from 
the GR11 sample, leaving a total of 60 galaxies. Two more of these galaxies are cut since 
they do not meet the assumed spatial resolution of our following simulations, leaving a total 
of 58 galaxies. Therefore, the data in our sample has an effective resolution of 0'.'08. The 
velocity dispersions range from 72 km s _1 to 335 km s _1 , while the black hole masses range 
from l.lxl0 6 M o to 5.6xl0 9 M P 
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2.2. Simulated Sample 



To construct a simulated sample, we start with the sa me method as iBatcheldorl ( 120101 ) . 
by querying the HyperLeda0 catalogue ( Paturel et al. 2003 ) for galaxies with measured values 
of the central velocity dispersion. The search is limited to galaxies within 100 Mpc, the 
same limit that we place on the GR11 sample. We use two different parameters to retrieve 
galaxies from HyperLeda: modz and modO. The first is redshift-dependent and yields galaxies 
further than ~ 8 Mpc away, while the second is redshift-independent and fills in the sample 
for distances down to ~ 0.7 Mpc. Combining these two samples gives complete coverage 
from very small distances all the way out to 100 Mpc. The redshift-independent distance 
measurements from modO only captures galaxies closer than ~ 50 Mpc, while modz fills in the 
sample with galaxies all the way up to 100 Mpc. Figured] shows the distribution of distances 
in the sample. The complete base sample contains 2, 870 galaxies from 0.7 — 100 Mpc, with 
60 < a < 400 km s _1 , the range of interest from the observed GR11 sample. This sample is 
referred to as HL, and all M bh are simulated as described in Section 3. 




20 40 60 80 

Distonce [Mpc] 



Fig. 1. — The distribution of distances in the HL sample. The solid line shows the distri- 
bution for galaxies selected with modO, and the dot-dashed line shows the distribution for 
galaxies selected with modz. 



3. Methods and Results 



We use Bayesian Monte-Carlo simulations to build models that can constrain and in- 
terpret the observed — a distribution. Based on Bayes' theorem, given the observed 
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distribution T> in the M^h — a plane and the combination of selection effects S, the proba- 
bility of a hypothesis P(T-L\T>,S) is proportional to 

P(H\V, S) oc P(V\n, S)P(U\S), (1) 

where P(1-L\S) is the prior probability of the model and P(T)[H, S) is the likelihood function. 

To construct intrinsic models Ti, we use Monte-Carlo simulations assuming three dif- 
ferent models for the black hole mass distributions: (i) an Mbh — a relation characterized 
by three parameters, a slope (3, an intercept a, and an intrinsic scatter e, (ii) a power law 
distribution P(M) oc M _r , and (iii) a power law distribution with an exponential cutoff at 
a fixed mass of 10 9 M Q . These models can be divided into two classes of model — a 
distributions. The first class, containing Model (i), is a simulation where the SMBH mass 
for each galaxy is calculated from the HL sample using an intrinsic Mbh — a relation, and 
therefore Mbh is dependent on a. The second class contains two simulations, Models (ii) and 
(iii), both of which generate M bh independently of a. The first class of models presupposes 
an intrinsic Mbh — a relation, and is the focus of this paper. The second class of models is 
used to explore whether selection biases alone are enough to recreate the observed Mbh — a 
relation, and to serve as a comparison for the first class, where Mbh is a variable that is 
dependent on a. By comparing the outcomes of the two classes, we can comment on which 
outcome is more likely, an intrinsic relation, or a relation that arises only because of selec- 
tion effects. In all cases, we use the HL sample with 2,870 galaxies for our a, and simulate 
SMBH masses for each galaxy. For each distribution, we focus on the mass range of 1O 6 M 
to 10 10 M Q . 

Once the intrinsic models are constructed, we apply selection effects, or selection func- 
tions S. The first selection function is the familiar sphere of influence cut («S rj ), which selects 
galaxies where the ratio of the sphere of influence to t he resolution is > 1. We u se a resolu 



tionof 0'.'08, which is close to the resolution of O'.'l that lGultekin et al.l (120091 ) and lBatcheldor 



(120 101 ) use. This provides a slightly looser constraint on the lower bound, but is still enough 
to identify implications from this selection effect. The second selection function applied, S a , 
is dependent on the velocity dispersion distribution. Ideally, the velocity dispersion distribu- 
tion of the observed sample should follow the general distribution for all galaxies to reduce 
the selection effect, especially if the underlying relation is non-linear or distorted. Because 
the general Mbh — & relation is measured for a combination of different type of galaxies, we 
use the velocity dispersion distribution from HL, a large sample of galaxies, as a proxy for 
our comparison. Figure [2] shows the velocity dispersion distributions in the H L and GR11 



sampl es, as well as the velocity dispersion function for early-type galaxies from ISheth et al. 



(120031 ). We find that in the observed sample (GR11) there are relatively more black mass 



measurements at the high velocity dispersion regime, and less in the low velocity dispersion 
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regime. A Kolmogorov-Smirnov (K-S) test comparing the GR11 and HL samples returns 
a p-value of 0.008, indicating that the samples are consistent only within 3a of each other. 
Comparing the HL sample and the velocity dispersion function for early-type galaxies, the 
velocity dispersion for early-type galaxies is narrower than the HL sample. Therefore, the 
velocity dispersion distribution in the observed sample for measuring the Mj,h — a relation 
does not follow that of a large sample of galaxies, and we include this selection effect in 
our simulation to test if this difference can cause significant changes in our final parameter 
estimations for intrinsic models. To match the distribution for the GR11 sample (applying 
the S a selection), we divide the a-space into five equal bins, and then S a cuts the number 
of galaxies in each bin from the model "H, such that the number ratios are consistent with 
those of the GR11 sample. The third selection function is a combination of both the sphere 
of influence cut and velocity dispersion sample bias, and is denoted as <S CT>rr 




Fig. 2. — Histograms showing the normalized distributions of velocity dispersion f or the 
HyperLeda (solid sha ded green area ), early type galaxy function from ISheth et al.l ( 12003 



hashed red area) and I Graham et al.l (120111 . solid black line) samples. The HL and GR11 
samples are only consistent within the 3a limit. 



After simulating a model and applying selection effects, we obtain a large number of 
scattered data pairs, and we compare the observed Mbh — a data pairs with the simulated 
data pairs to calculate the likelihood function P(T>\H, S). First, we use the two-dimensional 
K-S test with and without rotations to calculate the likelihood function. However, we find 
the K-S test is too general for this specific case, where the observed data fall in a linear 
relation, and the resulting constraints from our simulations are not tight. Since the sphere 
of influence selection causes the resultant simulated data points to always look linear, we use 
a linear model parameterized by a slope, an intercept, and an intrinsic scatter (/?', a', e') to 
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model the simulated data. Essentially, we assume P(D\H,S) = P(T>\/3',a',e') and evaluate 
P(T>\/3', a', e') using the standard y 2 likelihood including intrinsic scatter. Comparing the 
two methods, we find that the final expected values from posterior probability distributions 
are quite consistent; however, the uncertainties using the K-S method are always larger. 
Therefore, we subsequently only discuss results from the second approach. 

Since there is no prior knowledge on the model parameters given the selection effects, the 
prior P(H\S) is essentially proportional to the distributions of the model parameters, e.g., 
(3 in the — a relation, where we assume uniform distributions for the model parameters. 
To explore the three dimensional parameter space with the — a simulations, we use 
ranges of 7.8 < a < 9 for the intercept, 1 < < 10 for the slope, and 0.1 < e < 1 for the 
scatter. All of these ranges use a step size of 0.01, and cover the expected area of interest. 
For power law and exponential cutoff distributions, we assume a power law slope range of 
1.0 < T < 1.6. For model parameters outside of these ranges, the likelihood function is close 
to zero and will not contribute to the final parameter estimation. Finally, we measure the 
posterior probability distributions, P(H\U,S), using Equation 1, and estimate the expected 
model parameters as X = Pi^i/ X] P% anc ^ ^ a uncertainties from P("H|D, S). 

i i 

3.1. Model (i): Assuming an Intrinsic — a Relation 

We probe the entire three dimensional parameter space by letting a, (3, and e vary, 
calculate the x 2 likelihood of each individual combination of parameters, and construct 
the posterior probability distributions. We compute the expected value for the intercept, 
slope, and scatter X = ^PjXj/^Pj from the 3D posterior probability distribution. We 

i i 

marginalize the other two parameters for a particular parameter, X, by summing all the 
probabilities for a specific value of X (Figure [3]). Figure [3] shows the marginalized posterior 
probability distributions for each parameter for each application of selection effects, and 
where the expected values for each parameter of the simulation are marked. Integrating 
under each curve gives the 68% limits on the expected values. 
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7.9 8.0 8.1 8.2 8.3 3 4 5 6 7 8 0.2 0-3 0.4 0.5 0.6 0.7 

intercept, a slope, (S scatter, e 



Fig. 3. — Marginalized posterior probability distributions for the three parameters that 
describe the M^h — <J relation, a, (3, and e. The left column shows the expected values for the 
intercept a, the middle column for the slope j3, and the right column for the scatter e. The 
top row is the unfiltered simulation, the second row shows the simulation with the velocity 
dispersion selection bias («S CT ), the third row the sphere of influence selection effect (S n ), and 
the bottom row with both selection effects (S ajr .) applied. The green, solid vertical line is 
the expected value, and the 68% limits are depicted by the blue, dashed vertical lines. The 
summed probabilities have been normalized to have a peak value of unity. 
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Table 1. Best-fit Parameters of the — ex Relation 







Literature 






A/3 correction 


Reference 


oc 


P 


€ 


Sample size 


a n 

Pmtr ^measured 


Gebhardt et al. (2000) 


8.08 ±0.1 


3.75 ± 0.3 


0.30 


26 




Ferrarese & Merritt (2000) a 


8.14± 1.8 


4.80 ±0.54 




12 




Ferrarese & Merritt (2000) b 


8.77 ± 1.4 


5.81 ±0.43 




29 




Tremaine et al. (2002) 


8.13 ±0.06 


4.02 ±0.32 


0.25-0.30 


31 




Ferrarese & Ford (2005) 


8.22 ±0.06 


4.86 ±0.43 


0.34 


25 




Hu (20081 


8.18 ±0.06 


4.57 ±0.37 


0.42 


18 




Graham (2008) 


8.13 ±0.06 


5.22 ±0.40 


0.33 


50 




Giiltekin et al. (2009) 


8.12 ±0.08 


4.24 ±0.41 


0.44 


49(±18 upper limits) 


0.25 


Graham et al. (2011) 


8.13 ±0.05 


5.13 ±0.34 


0.43 


64 


0.76 


Beifiori et al. (2012) c 


8.19 ±0.07 


4.17 ±0.32 


0.41 


49 


0.19 


McConnell et al. (2011) 


8.29 ±0.06 


5.12 ±0.36 


0.43 


65 


0.50 


This paper 


o n7 +0.08 


5.28i^ 


n 45+0 08 


58 





a ri required to be resolved. 
"n not required to be resolved 

c Sample B of the paper, which uses secure estimates of Mbh rather than upper limits. 

Note. — Measured slope, /3, and intercept, a, for the — a relation in the form Io^M^/Mq) = a + 

^log(o-/200km s" 1 ) 




a. [km/s] „. [km/s] a. [km/ 5 ] 



Fig. 4. — Simulation results plotted with both selection effects applied. The center panel 
shows the simulation for the expected value of (3 = 5.28, while the left and right show cases 
of simulations (/3 = 4.74 and j3 = 6.15) corresponding to the 68% limits of the expected 
value. The black circles are from our simulations, while the red squares correspond to points 
in the GR11 sample. Information on the method used (triangles for masers, squares for gas 
or stellar kinematics) and resolution of the sphere of influence (filled symbols for resolved, 
open for unresolved) comes from G09 and FF05. 
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The expected values for the parameters are listed in Table [2j It is important to note that 
these are the expected values for the intrinsic parameters from posterior probability distri- 
butions, i.e., what is used to generate the simulation that populates the M bh — a space. When 
selection effects are applied, the measured parameters from a direct fitting method without 
modeling the selection effects can be different from intrinsic parameters. For example, for 
a given intrinsic slope, the directly measured slopes are lower than the intrinsic slope for 
j3 > 4. This can be seen in Figure [5j where the directly measured value f3' is plotted against 
the intrinsic slope /3 that generated the intrinsic model. The selection function including 
both selection effects has the largest influence on the intrinsic model. The selection effects, 
therefore, cause a meas urement of slope that is significantly lower than the intrinsic slope 



for > 4. For example, iFerrarese fc Fordl (]2005[ ) measure the slope to be 4.86 ± 0.43 when 



applying the sphere of influence resolution criteria, which our simulations show would have 



an int rinsic slope of /? ~ 6. This is close to the slope of 5.81 ± 0.43 that IFerrarese fc Ford 



(120051 1 measure without the sphere of influence selection effect, supporting the results of our 



simulations. 

To illustrate the correlation between the model parameters, for each 2D parameter pair 
we marginalize the third dimension and plot the confidence level contours in Figure [6J The 
confidence levels are found by integrating along isometric contours under the 2D surface. 
For example, the top left panel of Figure [6] shows a vs. with the third dimension, e, 
marginalized, for the base simulation without any selection functions applied. The expected 
values, a, (3 are marked by a star and the peak position of the distribution is marked with 
a circle. We then integrate under the surface from the peak position until reaching 68% of 
the total volume under the surface, providing the la contour. The contours become larger 
as selection effects are applied, increasing the amount of error that is introduced into the 
measurement of the — a relation. 

Table 2. — o Expected Values 





a 


B 


c 


H 

H,S a 


8 09+ 05 
o -,0+0.04 
S ' 18 -0.05 

8.071°;?* 


4 52+ ' 34 
-0.30 

A 79+0-33 

0.32 

5-42l°;S 
5.28l°j* 


34+ - 04 
u -° -0.03 

39+ 05 
u - oa -0.03 

u '^-0.07 

u - 4o -0.07 



Note. — Expected values from the 3D posterior 
probability space of the Mj,^ — a relation using 
Bayesian Monte-Carlo simulations. 
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3.2. Models (ii) and (iii): No M^ — a Dependency 

It is important to test whether the upper bound of the M^—a relation can be reproduced 
with a more general model between SMBH mass and bulge velocity dispersion, given the 
presence of the sphere of influence selection effect. Models (ii) and (iii) assume that the 
velocity dispersions and SMBH masses are independent variables. We draw the velocity 
dispersions from HyperLeda, and randomly assign them to simulated SMBH masses. We 
test two different SMBH distributions. The value of these tests lie not only in seeing if we 
can reproduce the upper bound of the M^h — o~ relation with SMBH mass distributions that 
are independent of a, but also as a comparison with Model (i), where an intrinsic relation is 
assumed. 

3.2.1. Model (ii): Power Law 

The first SMBH mass distribution is described by a power law distribution, P(M) oc 
M~ r . The power law index ranges from 1 < T < 1.6. We compare the probability function 
of this model with the peak probability of the linear — a model from our previous 
simulation, and plot the ratios in the top left plot of Figure [7J This ratio is shown for 
every value of T within our simulated range, for the base model without applying selection 
effects and the models with selection functions applied. The ratios, therefore, represent 
the likelihood of the hypothesis that velocity dispersion and SMBH mass are independent 
variables compared with the hypothesis that SMBH mass depends on velocity dispersion. 
It is immediately obvious from the top left panel of Figure [7J that the selection functions 
drastically increase the probabilities. There is a significant difference between selection 
functions that use the sphere of influence selection effect and those that do not; we plot 
them in individual panels since those without the sphere of influence selection effect applied 
are several orders of magnitude smaller. We further examine the simulation with selection 
function S a>ri applied, but the simulations are never able to reach the la point needed in 
order to be consistent with the GR11 sample (depicted in the first row of Figure UJ along 
with the expected value of the power law slope, f = 1. 45^05). The bottom panel of Figure [7J 
plots the simulation for T = 1.45 along with the GR11 sample. There are too many outliers 
in the simulation above the relation, and too few points at the low velocity dispersion end 
of the simulation. This does not accurately capture the behavior of the observed — a 
relation, and can therefore rule out the model of a power law mass distribution coupled with 
selection effects as an explanation of the empirical — a relation. 
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Fig. 7. — For power law (left) and exponential cutoff (right) mass distributions: the top 
panels show the x 2 -likelihood ratio of the power law or exponential cutoff distribution to the 
maximum x 2 -likelihood value for the linear — o relation simulation, plotted as a dashed 
green line. The top panels have the sphere of influence selection effects applied. The bottom 
plots show the GR11 sample overlaid with the %,S ar . simulation for the expected values 
f = 1.45 (power law) and f = 1.30 (power law with exponential cutoff). 



3.2.2. Model (in): Exponential Cutoff 

The last case considered is a power law distribution with an exponential cutoff. This 
is the same as Model (ii) but with an exponential cutoff set at 1O 9 M . The cutoff shifts 
the expected value to T — 1.30loi4 f° r the simulation with S air . applied. The right side 
of Figure [7] shows the results for the expected value. This model is also not able to reach 
a point where it is consistent with the GR11 sample. The clustering of simulated galaxies 
is too heavy on the low-mass side of the relation, and while the scatter is tighter than in 
the power law model, there are still many outliers in this area, and not enough in the high 
mass end of the distribution. Therefore, we can also rule out this SMBH mass distribution 
combined with selection effects as an explanation for the — a relation. 
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4. Discussion and Conclusions 

Using a Bayesian Monte-Carlo analysis, we constrain the parameters of the Mbh — <J 
relation as a = 8.07±g;?g, /3 = 5.28±o:ii, and e = QA5 -om- The main results of this paper 
are: (1) the directly measured slope of an Mbh — o relation is shallower than the intrinsic 
slope because of selection effects in and a, and (2) the uncertainties on the slope of the 
relation are likely underestimated in previous measurements. We additionally conclude that 
the Mbh — cr relation cannot be successfully reproduced with either power law or exponential 
cutoff mass distributions coupled with r, and a selection effects, and that it is statistically 
much more likely that the relation is intrinsic rather than observed. 

Comparing the two selection effects shows clearly that the sphere of influence selection 
effect dominates, which will decrease the slope if the intrinsic slope is steep, e.g., (3 ~ 10, 
or increase the slope if the intrinsic slope is shallow, e.g., /3 ~ 1, with the crossing point at 
/3 ~ 4 (Figure E]). The rj selection resembles a parabolic curve in the logarithmic Mj,h — & 
plane, which can be approximated with a linear cut of /? ~ 3.5. For a steep intrinsic — o 
relation, more objects with lower masses will be cut if one applies the r» resolution selection 
criteria, resulting a shallower relation. The situation will reverse for very shallow intrinsic 
slopes (e.g., (3 ~ 1), where the simulated slope will be steeper than the intrinsic input 
slope after applying selection effects. However, since the likelihood function for these cases 
approaches zero, this effect is not important. For the power law and exponential cutoff mass 
distributions, our models assume more galaxies with smaller Mbh given the range of T used 
in the simulations. However, in the high a range, the selection will trim most of the low 
mass objects, and in the low a range, only a few objects have measurements in Mbh and it is 
more likely for the mass to be small based on our assumptions. Even though we are able to 
rule out these two scenarios, it is still clear that selection effects can significantly increase the 
chance that these intrinsic power law models would be consistent with the observed sample. 

The a selection effect can be important if the intrinsic relation between Mbh and a is 
non-linear as in our power law models. For linear models, the a selection effect can also 
be important if the intrinsic slope is extremely deep, f3 ~ 10, or shallow, /3 ~ 1, where the 
effect can modify the slope by A/3 ~ 1. Fortunately, within the range of intrinsic slopes of 
P = 4-6, the a selection effect is not important. 

The argum ent for using an r, select ion criteria is that black hole mass estimates can be 



unreliable (e.g., iFerrarese fc Fordl 120051 ) . Ensuring that the sphere of influence is resolved 
is thought to more accurately measure the Mbh — & relation. In this paper, we find that 
this will cause selection effects that impact estimating the best-fit parameters, especially the 
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slope f3. Based on our simulations, we find 



intr 



4.69 -f- 2.22 X pmeasuredi 



(2) 



where (3i ntr is the intrinsic slope and f3 measure d is the directly measured slope without modeling 
the selection effects. In the relevant parameter space {3 = 4-6, the r» sele ction effect decreases 
the m easured slope. This is consistent with some previous results (e.g.. iFerrarese fc Merritt 
20001 ) . The authors measured the M^ — a relation twice, once with secure measurements and 



the sphere of influence resolved, and once for their entire sample. They found that the slope 
decreases from 5 .81 =fc 0.43 to 4.80 ± 0.54 when including only galaxies where T{ was resolved. 
Giiltekin et al.l (120091 ) used Monte-Carlo simulations to show that an selection criteria 
can bias measurements of the zero point, slope, and scatter. They consider three different 
scenarios using synthetic data sets. In the first scenario, they require that be resolved based 
on the values of Mbh and a in the data set. For the next scenario, they use a to calculate 
the expected Mbh, find the radius of the sphere of influence using a and the expected mass, 
and then cut any data pairs where is not resolved. The last scenario uses an observed 
sample of a with simulated masses. The results here agree with their first scenario, as the 
slope decreases with selection effects. However, we find an opposite bias in the slope for the 
third scenario, which is closest to the simulations here as it uses observed a. The differences 



might stem from the fact that the lGiiltekin et al.l ( 120091 ) scenario uses observed distances and 
instrumental resolution, which give different values for r i; rather than assuming a standard 
resolution for all data points, as done here. In addition, iGiiltekin et al.l ( 120091 ) studies a 
sample yielding a low slope, which, based on our simulation, has a smaller correction on (3 
(Figure [5] and Equation 2). To test the robust ness of our conclusions, we run our Bayesian 



Monte-Carlo code with other recent samples (IGiiltekin et al.l l2009r iMcConnell et al.l 12011 



Beifiori et al.ll2012l ). We list our corrections on the slope, A/3, in Tabled) We find decreases 
of the measured slope after applying selection effects in all the samples, however, with larger 
corrections for larger measured slopes. This is consistent with the results of our simulation. 

There are a number of measurements of the My, — cr relation, but we do not find an 
analysis that fully models various selection effects. The fitting results among different groups 
are not consistent, due to the mass measurements, sam ple selections, or line ar regression 
fitting techniques (about A/3 ~ 0.2). Using the sample of iGraham et al.l ( 1201 ll ) but filtering 



out the two galaxies not satisfying our simulation resolution of 0'.'08, we find a best estimated 



slope of (3 = 5.28^0 55 by modeling the selection effects. Although this is i n agre ement 



with studies that measure higher values for the slo pe (e.g., IFerrarese fc Fordl ( 120051 ). (3 



Graham! (l2008h . j3 = 5.52 ± 0.40), this 



5.81 ± 0.43; IGraham et al.l (EOllt ). (3 = 5.13 ± 0.34; 
consistenc y needs to be fu r ther e valuated. For example, our analysis indicates that directly 
fitting the IGraham et al.l ( 120111 ) sample should yield a slope of (3 ~ 4.5 (Table 2), and 
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indeed, we find a slop e of ~ 4.6 by directly performing a linear regression on the volume- 
limited |Grjliairi_el_alJ ( 1201 ll ) sample with which we compare our simulations. Applying the 
selection is important to keep the mass measurement reliable; however, we also need to model 
this effect to restore the intrinsic slope of the relation. For groups measuring a shallower 
slope, ~ 4, the selection effect is small based on our simulation, and for measurements of 
steeper slopes ~ 5, the selection effect is larger. Therefore, it appears the difference in the 
measurements of slopes (~ 4 vs. ~ 5) among different groups arises mainly from the basic 
mass measurements. 

An Mbh — o relation has far reaching implications, as the slope of the relation can illu- 
minate important physical processes connecting a central black hole to the host galaxy. The 
results presented here suggest that the slope is likely higher than is measured d ue to ry, and 
a selection effects, and closer to the energy-driven winds scheme as predic ted by ISilk Sz Rees 
(M o c a 5 ; 119981 ) rather than the momentum-driven winds predicted by iFabianl (M oc a 4 ; 
19991 ). This feedback mechanism can be realized by quasar winds, especially in broad ab- 
sorption line quasars (BALQSOs). The efficiency depends linearly on the covering fraction 
of BALQSOs, and recent studies suggest that the intrinsic fraction of BALQSOs in quasars 



is at l east two times higher than those fractio ns measured in optical surveys (e.g., iDai et al. 



20081 : IShankar et al.ll2008al : I Allen et all 1201 lh . For low-ionization BALQSOs (LoBALs), in- 
cluding ones exhibiting iron absorption (FeLoBALs), the intrinsic fra ctions are estima ted to 
be 5-7 times larger than the values obtained from optical surveys (IDai et al.1 l2010bl ). Al- 
though the intrinsic fractions of LoBALs are still small, the column density in these wind 
can be orders of magnitude higher than other quasar winds. Using more robust measure- 
ments of the column densities of FeLoBALs, it is estima ted that FeLoBALs can provide 
kinetic feedback efficiency Eunl Lboi of a few percent (e.g., iMorabito et al.l 120111). consistent 



with the requirement for the feedback efficiency of the — a relation (e.g., ISilk &: Rees 



1998). Although th e energy feedback model can suffer from radiation loss (e.g., Kingj 



2003 



Silk fc Nusserll2010l ). it is possible that the feed back process is a multi-phase process, wher e 



the AGN feedback is just the first phase (e.g.. iHopkins fc Elvis! l2010l ; ISilk fc Nusserl 120101 ). 
Subsequent feedback from supernovae and stellar winds will provide additional energy. The 
slope of th e predicted Mu, — cr relation is still oc a 5 if the feedback is dominated by energy 
exchanges ( ISilk fc Nusserl 120101 ). 



An alternative model that is not studied here is the combination of an upp er envelope 
and s election effects causing us to see only the upper portion of the M hh — cr plane ( iBatcheldor 
20101 ). as physically interpreted by ?. However, iGultekin et al.l ( 120111 ) argue that these upper- 
limit models predict detection rates of black hole masses that are not observed. The lower 
portion of the M^h — & plane will be probed as technology and techniques advance enough 
to be able to resolve low mass black holes. Continuing to expand the observed samples 
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with reliable SMBH mass estimates is essential, but until we reach the point where selection 
effects are negligible, they must be well modeled when analyzing a sample. 

Here we have only tested a limited number of hypotheses against the current data, and 
therefore our conclusions are limited to those hypotheses. There are other factors that are not 
explored: for example, many studies differentiate galaxies by morphology when measuring 
the M bh — a relation (e.g., Graham et al. 2011; Hu 2008; Beifiori et al. 2012). While there are 
tantalizing hints that the — a relation may be different for different galaxy populations, 
the sample sizes are still too limited to be conclusive. This paper provides a formalism that 
is capable of testing a large number of hypotheses combined with selection effects, a useful 
springboard for future, deeper studies. 
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Fig. 5. — The directly measured slope, 0', vs. the intrinsic /3 of the simulation. The curve 
with no selection effects has a slope of unity, as expected. The selection effects cause the 
directly measured slope to be different than that of the intrinsic model. 
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Fig. 6. — The shading represents different confidence levels: the darkest shading shows 
the 68% contours, and increases to the 99.73% limit at the lightest shading. The highest 
probability in a simulation is marked with a circle, and the expected values are marked with 
a star. The top row displays the contours for the simulation without selection effects, the 
second row the velocity dispersion selection bias (S a ), the third row shows the simulation 
with the sphere of influence selection effect (S ri ), and the bottom row is with both selection 
effects ((S^rJ applied. From left to right, the columns show (3 vs. a, e vs. a, and e vs. /3 for 
the expected value of the third parameter. 



